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Abstract 

We elaborate on the interpretation of some mixed finite element spaces 
in terms of differential forms. First we develop a framework in which we 
j /^ ■ show how tools from algebraic topology can be applied to the study of 

their cohomological properties. The analysis applies in particular to cer- 
tain hp finite element spaces, extending results in trivial topology often 
referred to as the exact sequence property. Then we define regularization 
operators. Combined with the standard interpolators they enable us to 
prove discrete Poincare-Friedrichs inequalities and discrete Rellich com- 
pactness for finite element spaces of differential forms of arbitrary degree 
on compact manifolds of arbitrary dimension. 
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1 Introduction 



A close kinship between some mixed finite element spaces and some old con- 
structs of algebraic topology such as simplicial cochains was pointed out by 
Bossavit [S] and has received an ever increasing attention among numerical an- 
alysts, see in particular Arnold [2J. Early papers on mixed finite elements include 
Raviart-Thomas [33], Nedelec [2H][3I] and Brezzi-Douglas-Marini [IT]. Surveys 
on mixed finite elements can be found in Brezzi-Fortin [12] and Roberts- Thomas 
|35j . whereas for material on simplicial cochains we refer to Spanier [38 and 
Gelfand-Manin [3D] . This link has led to a reinterpretation of many results in 
numerical analysis, see for instance Hiptmair |23j and Boffi [7 , and has inspired 
the construction of new finite element spaces, see for instance Arnold- Winthcr 
[3] and Buffa-Christiansen [T4] . 

More specifically, some interpolation operators are known to provide a com- 
muting diagram, linking the De Rham complex of smooth differential forms to 
complexes of finite element spaces of piecewise polynomial forms satisfying com- 
patibility conditions on interfaces between the cells of the mesh. Based on such 
commuting diagrams it is often remarked that when the De Rham sequence is 
exact the finite element sequence is also exact. To carry out this reasoning one 
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usually extends the interpolators from smooth forms to spaces containing also 
finite element forms, which leads to some technical difficulties. For instance 
one often considers Sobolev spaces of differential forms. But continuity of in- 
terpolators on such spaces is based on the Sobolev injection theorems (see e.g. 
Adams-Fournier [1]), which in higher dimensions require W sp spaces with s and 
p large. Such norm estimates are also used when the commuting diagram prop- 
erty is used to prove convergence estimates such as discrete compactness in the 
sense of Kikuchi [25] see Boffi [B][7]. This property is crucial to the study of 
eigenvector approximations, see Boffi et al. [5] and Caorsi-Fernandes-Raffetto 
[15] . The lack of continuity of the standard interpolators leads to technical dif- 
ficulties even in two and three space dimensions. For many problems in physics 
it would be desirable to have estimates in four space dimensions and dimensions 
higher than four arise for instance in applications to finance. 

The argument using commuting diagrams to prove that the discrete complex 
has the same cohomology as the continuous one, is limited to trivial topologies 
(trivial De Rham cohomology groups, except perhaps the first or last one). 
But for many problems in electromagnetics the topology of the device is non 
trivial and this is an essential feature of the way it operates. For instance in 
a transformer the electric currents of the coils excite an electric current on an 
iron torus, mainly along a certain harmonic vector field. The space of harmonic 
vector fields on a compact manifold can be seen as a realization of a cohomology 
group of (classes of) 1-forms via a Riemannian metric - this is part of Hodge 
theory, see e.g. Taylor pT3] (Chapter 5). Thus for the numerical simulation 
of such devices it is a minimum requirement that the cohomological properties 
of finite clement spaces are correct. Devices with non-trivial topology are also 
the object of recent numerical studies, see in particular Rapetti-Dubois-Bossavit 
[32] . However these authors restrict attention to lowest order finite elements and 
domains in R 3 with connected boundaries. Partial results on the dimensions of 
various finite element cohomology groups have long been obtained using Euler- 
Poincare formulas, see for instance Nedelec [30] , 

That the De Rham sequence on general manifolds has the same cohomology 
as lowest order finite elements is in fact a deep theorem of De Rham, since the 
lowest order finite elements correspond to standard simplicial cochains. This by 
now rather well-known correspondence between finite elements and simplicial 
cochains was pointed out by Bossavit 8 following a remark by Kotiuga and is 
detailed for completeness in Proposition 12.31 of this paper. For a proof of De 
Rham's theorem we refer to Weil [40 . In this paper we provide a sufficient con- 
dition under which other finite element spaces also have cohomology naturally 
isomorphic to the De Rham cohomology. The analysis applies in particular to 
certain hp finite element spaces of differential forms, whose efficiency for sim- 
ulating electromagnetic problems has been demonstrated by Demkowicz et al., 
see the review [T5] , 

In fact we will work with arbitrary simplicial complexes which need not be 
manifolds, so we shall rather prove that other finite element spaces have coho- 
mology isomorphic to simplicial cohomology. This generalization of the kind 
of topological space on which we construct functional spaces, is crucial to the 
method of proof. Simplicial complexes which do not generate manifolds also 
appear in applications: for instance, returning to the example of electromag- 
netic phenomena, antennas are often constructed by adjoining more than two 
metal sheets (called screens in this setting) along edges. Such antennas are 



also successfully simulated in many industrial codes, though there is a lack of 
theory for the error analysis of these simulations. This paper can also be seen 
as a contribution to a construction of a framework where such problems can 
be addressed. For instance there could be interest in extending the results of 
Buffa - Christiansen [13] to the case of several interconnected screens. We there- 
fore provide a detailed description of piecewise smooth differentiable forms on 
arbitrary simplicial complexes. 

We notice furthermore that a proof of the fact that certain hp finite element 
spaces have the right cohomology can be found in Hiptmair |24] (Theorem 3.7 
p. 273). Apart from greater generality, the main advantage of our result lies 
perhaps in the simplicity of the hypothesis and the efficiency of the proof, taking 
full advantage of the language and tools of homological algebra. Indeed the 
present results turn out to be a rather straightforward application of the basic 
tools of homological algebra and, if nothing else, this part of the paper might 
serve as an illustration of how these tools can be applied in a finite element 
setting. Probably the results of SJ2] and 331 will not come as a surprise to the 
reader, but the framework is perhaps at variance with those usually adopted by 
both algebraic topologists and numerical analysts. 

That the dimension of cohomology groups of various finite element space is 
the "right" one, can be seen as an algebraic stability property. But the anal- 
ysis of approximation of eigenvalues, of say the curl curl operator relevant to 
electromagnetics, also requires metric properties such as the discrete compact- 
ness property already referred to. In order to obtain such estimates it has been 
argued that one should look for interpolation operators which are projectors 
commuting with the exterior derivative, with enough continuity properties in 
terms of Sobolev norms. As already mentioned the standard interpolators lack 
suitable continuity. On the other hand so-called Clement interpolation [17] is 
defined on rough functions but is not inserted in a suitable commuting dia- 
gram. Variants designed to remedy on this have been designed, for instance in 
Bernardi-Girault [I] and Girault-Scott ^T\ but were not found to be directly 
transposable to our problem (though some of the techniques of proof we use 
are inspired by these papers). In this paper we construct other regularization 
operators. Unfortunately they do not commute with the exterior derivative and 
they do not leave the finite elements spaces invariant. However we can control, 
through an additional parameter, the amount by which these two properties 
fail, in appropriate norms. Composing with the standard interpolators gives us 
a tool sufficiently powerful for the goals we set for ourselves concerning stability 
and compactness properties of Hodge decompositions. 



Organization. The paper is organized as follows: In 32] we provide the def- 
initions of simplicial complexes we will work with and introduce a space of 
piecewise smooth differential forms satisfying a compatibility condition along 
interfaces. For brevity the latter will be referred to as compatible forms. We 
define lowest order finite elements as particular finite dimensional spaces of com- 
patible forms and check that they correspond to simplicial cochains (Proposition 
12. 3|) . Then, in SJ2J we prove that the inclusion of lowest order finite elements in 
the spaces of compatible forms, induces isomorphisms in cohomology (Theorem 
13. 7p . We then consider complexes which are intermediate between lowest order 
finite elements and general smooth compatible forms. Under natural assump- 



tions involving interpolation operators, we prove that inclusions (and interpo- 
lators) induce isomorphisms in cohomology (Proposition 13.81) . We also prove 
that the standard high order hnitc-elements behave well under wedge products 
(Proposition 13. 9p . Next, in 2] we restrict attention to smooth compact mani- 
folds (without boundary). We construct a regularization operator and use it to 
prove a discrete version of the Poincare-Friedrich inequality (Proposition I4.9P 
and Rellich compactness (Corollary 14. 131) . 

2 Definitions 

2.1 Simplicial complexes 

A simplicial complex is a set T of finite non-empty sets with the property: 

VTeTVT'cT T' ^ =*> T' £ T. (1) 

All simplicial complexes considered in this paper are themselves finite, so we 
shall take the liberty of not repeating this additional requirement. A classical 
reference on simplicial complexes is Spanier [38 chapter 3. A more abstract 
approach is exposed in Gelfand-M anin [5D] . The non-empty elements of T will 
be called simplexes, and a simplex with k 4- 1 elements will be said to be fc- 
dimensional or to be a fc-simplex. The elements of a simplex are called vertices. 
For each integer k we denote by T k the subset of T consisting of fc-dimensional 
simplexes. The set of vertices of (of elements of) T, which is the union of T, 
will be frequently be identified with T°. If T is a simplex and T" a subset of T 
we say that T" is a face of T ; if X" is /-dimensional we also call T" an /-face of 
T. 

For integer k the following notation is useful: 

[fc] = {0,1, ••-,£}. (2) 

Given a fc-simplex T, on the set of bijections [k] — > T, the relation defined by: 

r ~ t' -<=>• r~ t' is an even permutation of [fc], (3) 

is an equivalence relation which has two equivalence classes when k > 1. An 
orientation of T is the choice of such an equivalence class. We suppose that for 
each fc-simplex T we have chosen a bijection: 

ct t : [k] -> T - ( 4 ) 

We will equip T with the orientation induced by ctt (its equivalence class) . 
An affine realization of a simplex T £ T k is an injection: 

PT : T -► V, (5) 

into an affine space V such that the range px (T) has a fc-dimensional affine span 
in V. We allow V to have dimension larger than k. The closed convex hull of 
Pt(T) in V will be denoted \T\ and also called an affine realization of T when 
no confusion can arise as a consequence. An affine realization of a simplicial 
complex T is the data consisting of an affine realization of each of its simplexes; 
it is denoted |T|. If T, T" are simplexes in T such that T" C T there is a unique 



affine map X" — > \T\ which coincides with p T p T , on pt'(T'); it is injective and 
we call it the canonical injection : 

i TT , : \T'\ -»■ \T\. (6) 

We notice that for any T" C T" C T we have : 

i rr = idiy I and iTT" — iTT' ° It't" • (7) 

In the following we consider simplicial complexes equipped with one (and only 
one) affine realization. 

We denote by f2 fc (T) the collection of all families {ut)tet such that for 
each T € T, ut is a smooth differential fc-form on \T\ (by this we mean a k- 
form having a smooth extension to the affine space generated by \T\), and such 
that for each T,T' 6 T with the property that T' C T, we have the following 
compatibility condition on the pull-backs: 

{itt')*ut — Ut 1 - (8) 

The elements of $7 fc (T) will be simply referred to as compatible k- forms (or 
just compatible forms). The exterior derivative d and wedge-product A, ap- 
plied component-wise to compatible forms yield compatible forms since these 
operations behave naturally under pull-backs, see e.g. Lang [26 chapter V. 

If u — {ut)teT £ Q k {T)i it will be convenient to call ut the restriction of 
u to \T\ and denote it by u\t- If T' C T is also a simplicial complex we define 
the restriction of u to T' by: 

u\ T > = (ut)tzT' 6 n k (V). (9) 

It is sometimes convenient to associate with a simplex, the simplicial complex 
consisting of all its subsets. For instance if T is a simplex embedded in an affine 
space, this enables us to speak of J7 fc (T). We notice that if u € £l k (T) in the 
above sense, then u is uniquely determined by u\t, and that conversely any 
fc-form on \T\ gives rise to an element of J7 fe (T). We denote by dT the simplicial 
complex consisting of all strict subsets (called proper faces) of T. 

The following definition will be useful when we construct regularizations. 
For a given TeTwe denote by st(T) the star of T, defined by: 

st(T) = {T'eT : TnTV^}. (10) 

2.2 Finite elements 

We now define lowest order finite elements as finite dimensional subspaces of 
fi fe (T). This is essentially a reformulation of Bossavit [S] and the expressions 
appearing here can at least be traced back to Weil [40] . Weil (following DeRham) 
considers a good cover of a manifold and a subordinated partition of unity 
consisting of smooth functions. A cover gives rise to a simplicial complex: the 
nerve of the cover. We, on the other hand, suppose that the simplicial complex is 
given and use it to obtain a partition of unity consisting of compatible piecewise 
smooth functions. 

For each vertex i, denote by Xi the element of f2°(T) such that for each 
simplex T having i as vertex A,|t is the barycentric coordinate map on \T\ 



relative to i, and Xi is on simplices not having i as vertex. In other words 
Xi is the unique piecewise afhnc function on |7~| such that Xi(j) = 5ij for all 
i, j G T°. The linear span of the family (Aj)j G 7-o is denoted A°(T). 

For each integer k > 1 and T G T k , denote by At the compatible fc-form: 

k 
X T = k\ ^(-l) i A <7T(i) dA CTT(0 ) A • • • (dA ffT(l) ) A • • • A dA CTT(fe) , (11) 

i=0 

where the symbol (-) A signifies omission. The linear span of the family {X^TeT 1 * 
is denoted A fc (T). It can alternatively be characterized as: 

k 

span{^(-l) i M,du A • • • (dui) A • • • A du fe : Vi e [fc] m G A°(T)}. (12) 

i=0 

For any k we denote by V[k] the group of permutations of [k], and by e : 
V[k] — > { — 1,1} the signature morphism. We notice that for any functions 
«0) • • • ! u k G £1°(T) we have: 
k 
fc!^(-l)*Ujdu A- •• (du l ) A • --Adufe = ^ e ( /r ) u r(o)du T (i)A- • -Adu T ( fe ). (13) 

j=0 r£P[i] 

We define the degree of freedom associated with a fc-simplex T as the linear 
form Ht on O fe (T) defined by : 



[i T : u^ j u\ T , (14) 

\T\ 



when \T\ is oriented by cr^. 



Proposition 2.1 For ant/ iwo k-simplexes T and T', we have ht'{Xt) = o~tt' , 
where the last symbol is the Kronecker delta. 

- Proof: First we remark that, since the family (A,),gT constitutes a partition 
of unity on \T\, we have: 

(X t )\t = k\ (dA CTT(1) A • • • A dA ffT(fe) )| T . (15) 

The integral of this fc-form on \T\ can be computed in coordinates and is 1 by 
the choice of orientation and normalization. 

Next we remark that if T' ^ T we can pick an i G T \ X". The function Aj 
is on \T'\ (as well as its exterior derivative), implying that Xt\t' = 0. This 
completes the proof. D 

Since this proposition implies linear independence of the family {Xx)TeT k 
we obtain: 

Corollary 2.2 The family {Xr)T^T k * s a basis for A k (T). 

For any (k + l)-simplcx T and any fc-face T' of T we define the incidence 
number e(T,T') G { — 1, 1} to be the sign of the permutation of [k + 1] defined 
as: 



^ cr T (0), (16) 

?; + l ^ o- T l {o- T >{i) + l) for i G [k]. (17) 

For all simplexes T, T' not covered by this definition we put e(T,T') = 0. 



Proposition 2.3 The exterior derivative d maps A k (T) into A k+1 (T) and the 
matrix of d : A k (T) — > A k+1 (T) in the bases (Ajv)T'eT fc ~ * {^T)reT k+1 has 
entry e(T,T') at the indices {T,T') G T k+1 X T k . 

- Proof: First we remark that: 

dA T < = (k + 1)! dA CTT , (0) A • • • A dA CTT , (fe) . (18) 

Next we make the following computation. Let (iti)jgi denote any family of 
functions in fi°(7") constituting a partition of unity of \T~\ indexed by a finite 
set / containing [k] as a subset. Choose a k < n. We have: 

duo A • • • A duk 

= 2, Ujduo A • • • A duk + 2, u iduo A • • • A duu 
j#[k] ie[fc] 

= >_, Ujdwo A • • • A d?ife + N^ "id-Mo A • • • (d(l — \, u j)) • " ' ' A dw^ 

M*] «6[fc] M*] 

= ^ ( Ujduo A ■ • • A du k + J^ (— l) J+1 <u,dUj A du A ■ • • (du t ) A ■ • • A dw/c J 

If we change the index tuple (J, 0, 1, • • • , k) to (0, 1, • • • , k + 1) and remark that 
for any family (ui)ie[k+i] °f elements of J7°(T) we have: 

u dui A • • • A dM fe+ i + 2^2 (-l) !+lM »+iduo A • • • (du i+ i) A • • • A du fe+ i 

i£[k] 

= Yl {-iyuidu A---(du l ) A ---Adu k+1 , 
ie[A+i] 

then, comparing the signs appearing here with the definition of incidence number 
gives: 

dA T - = J2 t(T,T')\ T . (19) 

Ter k + 1 
This completes the proof. □ 

Remark 2.1 Suppose that the set T° has been given a total ordering and that 
for each T G T k one chooses the bijection o~t '■ [k] — > T to be the one which is 
increasing. Pick T G T k+1 and a k-face T" of T . Let I £ [k + 1] denote the 
integer such that: 

{a T (l)}=T\T f . (20) 

Then: 

e(T,T') = (-l) 1 . (21) 

We thus have a diagram : 

-> a x (n -> tt 2 (r) -> ■•• 

(22) 
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where the horizontal arrows represent exterior derivatives and the vertical ar- 
rows are inclusion mappings. We refer the reader to Appendix \K\ for some no- 
tions of homological algebra which will be applied to the above diagram. Since 



dod = 0, each row is a complex. The vertical arrows of course commute with the 
exterior derivative, providing a morphism of complexes. The cohomology group 
R k n'(T) is defined as the quotient of the kernel of d : n k (T) -»• Vl k+1 (T) by 
the range of d : Cl k ~ l (T) — > il k (T). The cohomology of A* (T) is defined simi- 
larly, and one checks that morphisms of complexes induce morphisms between 
cohomology groups. 

Remark 2.2 By Provosition \2.3\ and Remark \2.1l it follows that the complex 
A* (T) is isomorphic to the simplicial cochain complex defined for instance in 
Gelfand-Manin [201 %I.^.l. Said differently, the simplicial cochain complex is 
the coordinate expression of A*(7~) in the bases (A^reT* when the simplexes 
are oriented in a certain way. 

So-called interpolators n^- : il k (T) — > A fc (T) are constructed by assigning 
to u £ VL k (T) the unique element v £ A fe (T) such that for each T 6 T k : 

Urn = /it"- (23) 

Evidently, these are projections. Moreover we have: 

Proposition 2.4 The interpolators IT^- commute with the exterior derivative, 
i.e. we have commuting diagrams: 

n k (T) d n k+1 {T) 

n k T | n^ +1 | (24) 

K k {T) d A fc+1 (T) 

- Proof: Application of Stokes theorem in a simplex T, taking into account 
relative orientations, gives: 

[ d\ T > - $>( T ' T ") / A ^'> ( 25 ) 

= e{T,T'). (26) 

By Proposition 12. 31 ( specifically equation (fi~9|) ) this gives the desired result. D 

3 Cohomology 

We will now show that the inclusion mappings (vertical arrows) in (|2"2")l induce 
isomorphisms between cohomology groups. Then we will examine intermediate 
complexes appearing in finite element theory. 

3.1 An isomorphism in cohomology 

We examine first the properties of a single simplex. 

First we recall the construction of homotopy operators. Their efficiency in 
the context of finite elements was noticed by Hiptmair |23] . 

Let T be a p-simplex equipped with an affine realization pt : T — > V. Let 
xq be one of its vertices and define functions F t : V — > V by: 

F t (x) = x + t(x-xo), (27) 



Define also a vector field lonl^ by: 

X(x)=x-x . (28) 

Introduce finally A, a map taking fc-forms on \T\ to k — 1-forms on \T\, defined 
by (here v\X is the contraction of the form v by the vectorficld X): 

Av= I F*(v[X)dt. (29) 

Jo 

Then it is standard procedure to remark that for any function u on \T\: 

u — u(xq) = Adu, (30) 

and if u is a fc-form on \T\ with k > 1: 

u = Adu + dAu. (31) 

As a consequence we obtain: 

Proposition 3.1 We have dimH ft*(T) = 1 and for k>l, R k n'(T) = 0. 

As noticed by Hiptmair (even for higher order polynomials), by the special 
choice of homotopy F t , A also maps A k (T) into A k ~ 1 (T), and it follows that 
dimH°A*(T) = 1 and for k > 1, H°A*(T) = 0. For the purposes of future 
reference we therefore state the trivial consequence: 

Proposition 3.2 The inclusions A fe (T) — » CI (T) induce isomorphisms 

H fe A # (T) ^H fe fi*(T). 

Proof: Only the case fc = needs to be considered, and it is straightforward. 
□ 

We will also need a similar result for the boundaries of simplexes. First we 
state the following: 

Proposition 3.3 The restriction map fl k (T) — > Vt k {dT) is onto, i.e. compati- 
ble piece-wise smooth k-forms on boundaries of simplexes can be extended to the 
interior. 

- Proof: For i G [k] put Xi = pto~t{i)- In the following T will often be 
identified with the set {xq, • • • , Xk}, enabling us for instance to call the points 
Xi the vertices of T. Let U be the affine subspace of V generated by \T\ and 
Ui be the complement in U of the afhne span of the face of T opposite to x^. 
Let (fa) denote a partition of unity of a neighborhood of \T\ in U, subordinated 
to the open covering (Ui). Thus cf>i is 1 on a neighborhood of Xi and on a 
neighborhood of the affine realization of the face of T opposite to Xi . 

Pick an Z-face T of T and denote by T" the opposite face T\T'. Let 
Xi , ■ ■ ■ ,Xi, denote the vertices of T". In this setting we (temporarily) denote 
the barycentric coordinate on T associated with the vertex x% by Xi, and define 
Ft' to be the map \T\ \ \T"\ —> \T'\ defined in barycentric coordinates by: 

(A ,--- ,\ p )^(X l0 ,--- ,Ai,)/(Ai +-" + Ai,). (32) 

We also define <j>T' to be the function: 

<h' =4>i -\ 1- 4>%i ■ (33) 



Remark that <f>T' is 1 on a neighborhood of \T'\ and on a neighborhood of 
\T"\. 

Pick u G n k (dT). If V is a face of T and the restriction of u to <9T' is 
then the following fc-form on \T\ : 

E T 'U = 4>t'F%,{u\t>), (34) 

satisfies : 

(i TT ,)*E T ,u = u\ T >, (35) 

and for any other face T" of T with the same dimension as T" we have (since 
u\st' = 0): 

(i T T>>)*E T ,u = 0. (36) 

It follows that if u G fl (dT) is on all geometric realizations of faces of dimen- 
sion < I, then the restriction to dT of the fc-form on \T\: 

J2 E T ,u, (37) 

T ' e r' +1 

coincides with u on all those with dimension < I + 1. 
Now define an extension of u recursively by: 

E u= Y^ E T'U, (38) 

T'eT" 

and for < I < k - 2: 

E l+l u= Y. E TI {u-{E u+--- + E lU )\ dT ). (39) 

Then by the above considerations Equ + ■ ■ ■ + Ek~iu is indeed an extension 
ofuto|T|. D 

Remark 3.1 The similar result for A* follows easily from the fact that the 
barycentric coordinate maps Xi can be extend naturally from the boundary to the 
simplex (yielding barycentric coordinate maps). 

The space of fc-forms on T whose restriction to the boundary dT is is denoted 
Q k (T : dT). The exterior derivative commutes with restrictions yielding maps: 

d : Q k (T, dT) -^Q k+1 (T, dT) . (40) 

Thus we obtain a subcomplex of fl m (T) denoted fl'(T,dT). The next proposi- 
tion is about its cohomology groups (called relative cohomology groups). 

Proposition 3.4 For any k < dimT, R k n m (T, dT) = i.e. for any u £ n k (T) 
such that u\gx — and du = 0, there is v G £l k ~ 1 (T) such that v\qt = and 
dv = u. 

For k = dimT we have dimH k H,'(T, dT) — 1 and, for any u € il k (T) such 
that u\qt = and du = 0, there is v G Vt k ~ l {T) such that v\qt — and dv = u 
if and only ifj / u = 0. 



1 Integration being taken with respect to any chosen orientation of T. 

10 



The following induction argument is a variant of the one used to prove a sim- 
ilar result for polynomial forms in |23j (Lemma 17, which however does not 
distinguish the case k = dimT). 

- Proof: We proceed by induction on the dimension of T. For dimT = 
the result is clear. Suppose now that the proposition is true for simplexes of 
dimension < j for some integer j, and consider a simplex T of dimension j + 1. 

First we suppose that k < dimT. Pick u e il k (T) such that u\qt = and 
du = 0. We use the homotopy operators constructed above, relative to a vertex 
x Q of T. Let T' denote the face of T opposite to x - Put v = Au <E fi fc_1 (T). 
Then we have dv — u and for any proper face T" of T - except perhaps T' - 
v\t" = (since T" is in a proper face of T containing xo). We correct for this 
eventuality using the induction hypothesis: 

Put v' = v\t' ■ Then v' is a (k — l)-form on the (dimT — l)-simplex T', such 
that v'\dT' = and dv' = u\ T > = . We can therefore pick w' e £l fe ~ 2 (T') 
such that w'\dT' = and du/ = v' . Then extension of w 1 by to dT yields an 
element of fl k ~ 2 (dT) still denoted w', and we can extend this to an element w 
of fl k ~ 2 (T) by Proposition 13.31 Then v — dw satisfies d(v — dw) = u and by 
construction (v — dw)\gr = 0. 

Now we treat the case k = dimT. If du = u and v\qt = 0, the integral of u 
is necessarily by Stokes formula: 



u= dv= v\ dT = 0. (41) 

T JT JdT 

In particular any constant non-zero element J7 fc (T) provides a non-zero relative 
cohomology class. On the other hand if du = and J u — we can use a 
construction similar to the preceding case, noticing that, with (as above) u = Au 
and v' = v\t'- 



v' = / v\ 9T = dv= u = 0. (42) 

IT' JdT JT JT 

This completes the induction step and hence the proof. □ 

Proposition 3.5 The inclusions A k (T,dT) — > fl k (T,dT) induce isomorphisms 
H fe A*(T,aT) -> R k Q*(T,dT). 

- Proof: Pick u G A fc (T, <9T). We remark the following: 

- If k < dimT then u\qt = implies u = 0. 

- If k = dimT then J u — implies u = 0. 

From these two remarks and the preceding proposition the result follows. □ 

Proposition 3.6 The inclusions A k (dT) — } tt k (dT) induce isomorphisms 
R k A*(dT)^}i k n'{dT). 

- Proof: For each k we have a short exact sequences of the type: 

-» fl k (T, dT) -> fi fc (T) -^ fi fe (9T) -^ 0. (43) 

They provide a long exact sequence in cohomology (see appendix [AJ . We have 
a similar result for A*. Thus, for each k, we have a commutative diagram: 

R k n'(T,dT) ^H fe fr(T)-> R k n*(dT) -> E k+1 CT(T,dT) -> H fe+1 0*(T) 

H fc A*(T,9T) ^H fc A'(T)^ H fc A*(9T) -)• H fe+1 A'(T,aT) -> H fc+1 A'(T) 

(44) 
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where each row is an exact sequence, part of the long exact sequence in coho- 
mology obtained from the short ones, and the vertical arrows are maps induced 
by the inclusions of complexes A* —►£]!". 

In diagram (PHI) the two first and the two last vertical arrows are isomor- 
phisms. It then follows by the five lemma (see Gelfand-Manin [20] p. 120) that 
the middle arrow is an isomorphism. □ 

The preceding results on simplexes can be extended to simplicial complexes 
as follows: 

Theorem 3.7 The inclusion mappings A fc (T) — > £l k (7~) induce isomorphisms 
in cohomology H fc A # (T) ->■ H fe Q*(T). 

- Proof: The simplicial complex T can be constructed by starting with a 0- 
simplex and adjoining at each step a top-dimensional simplex whose boundary 
is already included. The case of a complex consisting only of a 0-simplex has 
already been treated (if need be) so we consider now a simplicial complex T' 
for which the theorem is true and a simplex T £ T' such that dT C T' and the 
dimension of T is larger than, or equal to, the largest dimension of the elements 
of T' '. Let T denote the simplicial complex T' U {T}. 

We consider the short exact sequences of Mayer- Vietoris type (see [38] chap- 
ter 4, §6): 

o -»• n k (T) -> n fc (T) x n k (T) -*■ n k (dT) -»■ o, (45) 

where the second arrow is the injection: 

u^(u\ T ',u\t), (46) 

and the third arrow is the surjection (surjectivity follows from Proposition 13 .31) : 

(u,v) *-tu\aT- v\dr- (47) 

These provide a long exact sequence in cohomology. 

We also have similar short exact sequences for A* providing a long exact se- 
quence in cohomology. The inclusion maps A* — >• CI' then provide a commuting 
diagram relating these two long exact sequences, much as in Proposition l3.6l 

H fc fi*(f ) x R k n'{T) -»■ K k n m (dT) -► H fc+1 0'(T) -> R k+1 n'(V) x R k+1 n m (T) -> H fc+1 n*(9T) 

H fc A'(T') x H fc A'(T) -^ H fe A'(c)T) -*■ H fe+1 A'(r) -»■ H fe+1 A'(r') x H fe+1 A'(r) -)• H fe+1 A'(aT) 

(48) 
Once again the five lemma then gives the desired result. □ 

3.2 On high order finite element spaces 

Consider the following situation: suppose we are given spaces T k (T) such that 
A k (T) C T k {T) C n k (T) and d maps T k (T) into T k+1 (T). We suppose 
furthermore that the spaces r fe (7~) are equipped with interpolation operators 
n r : fi fe (T) -> T fc (T) i.e. projectors with range T k (T), and that the inter- 
polators commute with the exterior derivative. We also suppose that we have 
equipped A k (T) with interpolation operators II a : Vt k (T) — > A fc (T) commuting 
with the exterior derivative. 

12 



Proposition 3.8 Suppose that IIa ° Ilr = IIa. Then the inclusions A fe (T) — > 
T k (T) induce isomorphisms in cohomology. 

- Proof: Since IIa is a left inverse of the inclusion map A* (T) — > &'(T) and 
the latter induces isomorphisms in cohomology, IIa also induces isomorphisms 
in cohomology. Thus the kernel complex (see appendix |X| of IIa is exact. 

Consider the morphism of complexes T* (T) — > A* (T) obtained by restricting 
IIa ! we shall prove that its kernel complex is exact, from which it follows that 
it induces isomorphisms in cohomology. 

Suppose u 6 T k (T) satisfies Titvu = and du = 0. We can pick a v <E 
fi 1 (T) such that H\v = and dv = u. Then Ilrw is an element of T k ~ 1 (T) 
such that riAiIrw = IIa.v = and dllrw = Ilru = u. 

This completes the proof. □ 

For lowest order finite clement spaces we have already constructed interpo- 
lators determined by equation (|23|) . For high or variable order finite elements 
the construction of suitable interpolators is more technical. Usually they are 
constructed by imposing conditions similar to (|23[) . for a given choice of sim- 
plexes T with dim T > k and where v and u are perhaps differentiated and then 
wedged with an appropriate space of polynomial forms before being integrated; 
one says that one imposes high order moments to match. Interpolators of this 
type, commuting with the exterior derivative, were constructed in Demkowicz 
et al. [IS], see also Hiptmair [51] §3.5. In any case we remark that most often, 
to construct interpolators for high order finite elements, one imposes matching 
moments with respect to a family of linear forms on J7 fc (7~) containing (l^T)TeT k 
as a subfamily, and this guarantees that the hypothesis IIa ° Ilr = IIa holds. 

One way of constructing higher order finite element spaces is the following. 
Put X k = A k (T) for each k. Define A° to be the vectorspace generated by 
functions of the form u\U2---u n with Ui G A" for each i. Then define, for 
k > 1 and n > 2, X k to be the vectorspace generated by compatible k- forms of 
the form uv with u £ X^ l _ 1 and v £ X k . Our main result on these spaces is: 

Proposition 3.9 For each k,l and m,n the wedge of forms provides a map: 

A : X k n xX l n ^ I*t (49) 

- Proof: We prove that A : X k x X[ — > .AT* , from which the proposition 
follows immediately. 

We will need some notations: 

Let (v,i) be some family of functions indexed by consecutive integers. For 
any set of consecutive integers k < ■ ■ ■ < I we put: 

Su[ k ...i] — du k A • • • Adui. (50) 

This notation will also be used when one index is missing in the set {k, ■ ■ ■ , I}, 
e.g.: 

& u \k.:i-i] = du * A ' " ( du A ' ' ' A du/. (51) 

We also put: 

l 

i—k 

a notation which is extended straightforwardly to the case of one missing index. 
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We will prove, by induction on k, that: 

fc-i 



u[ -k-i] /\u [k ... k+ i] = (-l) k 1 ^( -1 )' u < u [o-i-fc+l]- ( 53 ) 

i=0 

It is evidently true for k = 1 and, if it is true for a given fc > 1, we can make 
the following computations. We remark that: 

u[-i---k-i] Au[k...k+i] = (u-i<$U[o-k-i] -£«-i Aup.-.fc-i]) AU[fc... fc+ q. (54) 
Concerning the first term on the right hand side, we see that: 

fc-i 
u_i<fa[ -fe-i] A u [k ...i] = (-I) fc w_i(w [0 ... fc+;] - ^ w^w [0 ...5... fc+ ;]). (55) 

i=0 

For the second term we remark that (by the induction hypothesis): 

8u-i A U[ ...fe_i] A U[ k ... k +i} (56) 

fc-i 

= (-f) fe - 1 ( 5 W _ 1 A^(-fru lW[0 ...,... fc+i] (57) 

4=0 
fc-1 

= (-l)*" 1 Z)(-l) < «i(-«[_ 1 ...i... fc+ i]+tt-i*tt|p...i... fc+ q) (58) 

i=0 

Now the last term in (|55|) cancels with the last term in (J58I) and we are left with: 

fc-i 
W[-i...fe-i] A u [fc ... fc+/ ] = (-l) fe (w-iu [0 ... fc+i] - 5Z(-l)*WiW[_i...j... fc+ j]). (59) 

j=0 

This completes the induction and hence the proof. □ 

4 Regular izat ion 

Unfortunately the interpolation operators such as ±1^- do not have sufficient 
continuity properties with respect to Sobolev norms, a problem which becomes 
more acute in higher dimensions. To remedy on this we construct regularization 
operators and show some applications to discrete Poincare-Friedrichs inequali- 
ties and discrete Rellich compactness. 

4.1 Smooth manifolds 

Let T be a simplicial complex and V an affine space. Suppose we have a map 
p : T° — > V and define for each TeTa map pt ■ T — > V by pr{i) = p{i)] 
it is an affine realization of T iff the affine span of pt(T) has dimension the 
dimension of T. We suppose now that this is the case for each T. Then we have 
an affine realization of T . We say that p is non-degenerate if in addition, for 
each T, V e T: 

|T|n|r'| = |rnr'|. (60) 

In these circumstance we also say that p determines an embedding of \T\ in V. 
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Consider a Riemaimian manifold M , which for simplicity we suppose to be 
smooth, compact and without boundary. The dimension of M is denoted n. 
The Riemaimian metric gives rise to natural scalar products of forms on M . 

We suppose that we have built a simplicial complex T, with an affine real- 
ization, and for each T € T an embedding Et : \T\ — > M. We require that for 
each TT'eTifT'cT we have: 



=T' 



3^ ° ITT' ■ (61) 



We suppose furthermore that for all T, T" £ T the following condition on over- 
laps holds: 

S r (|T|)nS 7 v(|T'|)=S r »(|Tnr / |). (62) 

Finally we suppose that M is generated by T, i.e.: 

M = U Ter HT(|T|). (63) 

In fact we consider a countable family {Th) of simplicial complexes, indexed 
by a parameter h representing the mesh width of Th- This parameter runs 
trough a set of positive reals accumulating only at 0, and we are interested in 
the limit h — > 0. For notational simplicity we suppose that these simplicial 
complexes are two by two disjoint (the only reason is to assure that a given 
simplex "remembers" which simplicial complex it belongs to). 

Given a complex T'(Th) of compatible differential forms, as in the preceding 
section, one can transport it to a complex of forms on M denoted r*(M), 
by pulling back by the maps S^ for T G Th- These transported forms are 
L 2 (M) (with respect to the Ricmannian metric) and their exterior derivative in 
the sense of distributions coincides with the transported exterior derivative we 
defined on the affine realization of Th, and is also L 2 (M). 

As is standard we denote by fi fe (M) the space of smooth fc-forms on M, 
whereas the space of piecewise smooth ones which are compatible with Th is 
denoted fi£(M). 

4.2 Construction of regularizations 

To fix notations we recall the construction of regularization by convolution. Let 
Lp denote a smooth non-negative function on 1" with support in the open unit 
ball, whose integral is 1. For <5 > define tp s by : 

i P s (x)^S^ n ip{S" 1 x). (64) 

Thus if has integral 1 and support in the open ball centered at the origin and 
with radius S. We let R denote convolution by cp, acting on differential forms 
on W l . That is: 

(R s u)(x)= [ip s (y)u(x-y)dy, (65) 



This regularization will be transferred to M in three steps. In accordance 
with widespread usage we will call the star of an element of 7^™ (the star was 
defined in equation (1101) ). a macroelement. We suppose that there is a finite 
family (<Si)ie/ of rt-dimcnsional macroelements embedded in M. n such that for 
each h and each T € Th we can choose a piecewise affine homeomorphism: 

$ T :\Si\->\8t(T)\. (66) 
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These will be called reference macroelements. For each i E I we denote by Tj 
the central simplex of Si so that Si — st(Tj). Let tpi denote a smooth non- 
negative function on W 1 with support in the interior of \Si\ yet equal to 1 on a 
neighborhood of \Ti\. 

First wc define regularization in the reference macroelements. Define a reg- 
ularization operator Rf on \Si\ by: 

Rfu=(l-<ipi)u + R 5 {ipiu). (67) 

This operator regularizes uniformly on \Ti\ and does not deregularize anywhere. 
We will always suppose that 8 is so small that the 8 neighborhood of supp tpi 
has adherence Ai in the interior of \S%\. Then Rfu — u has support in Ai for all 
u. 

Then we define regularization on the macroelements of Th- For any T 6 7^\ 
define Rj* by: 

R s T u = ($T 1 )*'^i (*t)*«- (68) 

This notation is somewhat abusive, but ^r{Ai) is in the interior of | st(T)| and 
($ T 1 )*i?f($ T )" r as an operator L 2 (|st(T)|) -> L 2 (|st(T)|) does not modify its 
argument outside &x(Ai); we can therefore extend it outside $t(^4*) by the 
identity operator. 

Finally we order T£ , so its elements can be written T± , T^ , ■ ■ ■ , TV and 
define regularization R s h on M by: 

R s h = R s TN ...R s T2 R s Ti . (69) 

(The operator R s h depends on the chosen order). 

4.3 Properties of the regularizations 

Approximation results. For the rest of this section we consider a fixed (ar- 
bitrary) k. Unless explicitly stated otherwise we deal here with differential forms 
of degree k. We put X h = T\{M). 

A key lemma is the following local stability result: 

Lemma 4.1 For each 8 there is a constant C such that for all h, all T G 7^" 
and all u 6 L 2 (st(T)) we have: 

\\n h R{ u \\v {T) < c|M| L2(st(T)) . (70) 

- Proof: If we pull back to the reference macroelement Si assigned to T by 
the map $t, and consider first the action of R s h we see that we have several 
contributions corresponding to the different terms in the product (|69p . We have 
the action of Rf , which smoothes out on a neighborhood of |Tj|. But we also have 
the action of R» coming from the reference macroelements Sj corresponding the 
neighboring simplexes of T in 1~£ ■ Since the transition maps between reference 
macroelements are continuous piecewise affine, these actions do not destroy the 
property of being piecewise smooth and compatible on a neighborhood of \Ti\. 
Moreover the number of different possibilities for different macroelements to thus 
act upon each other is finite. For each integer r we therefore have estimates of 
the form: 

\\R{u\\ C r m \) < C1| U || L2(8t(Ti)) (71) 
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Then we use the continuity of LT^ : C r (|Xi|) — > L 2 (|Ti|) and transport the 
estimate back to M. □ 

For each 8 this lemma immediately implies uniform boundedness of the maps 
U h R s h : L 2 (M) ->■ L 2 (M). But we actually have more: 

Theorem 4.2 For each 8 and for each u € L 2 (M) we have: 

lim|| U -n /l i?^|| L2(M )=0, (72) 

- Proof: Since the maps Ilfcj?^ : L 2 (M) — > L 2 (M) are uniformly bounded it is 
enough to prove the above estimate for a dense subset of L 2 (M). 

If u is a pullback from an open subset U C K™ to M of a differential form 
which is the product of an alternating multilinear map on W 1 by the char- 
acteristic function of a compact subset of U with nice enough boundary (say 
Lipschitz), then the norm convergence follows from the remark that one can 
ensure convergence in the reference macroelement for constant functions, and 
bound the contribution of the macroelements touching the boundary of suppu, 
by the preceding local lemma. □ 

When u is in a space compactly embedded in L 2 (M) we can obtain uniform 
convergence estimates, as follows. 

Proposition 4.3 For each Banach space X which can be considered as a sub- 
space of L i (M) such that the canonical injection X — > L 2 (M) is compact we 
have, for each 8: 

lim||7-n h ^|| x ^ Ij2(M) = 0. (73) 



- Proof: This follows from the preceding theorem, using Lemma IB. II in the 
appendices, letting (A n ) be the sequence (ILkR^) indexed by h for fixed 8, and 
K be the injection X — > L 2 . □ 

The Galerkin spaces spaces Xh are not invariant under H^R^ - but almost: 

Proposition 4.4 For each e there is 6' such that for each § < §' we have, for 
all h and all Uh S X^: 

\\u h - U h Riu h \\ L 2( M) < e||u/ l || L 2 (M) . (74) 

- Proof: For each reference macroelement the space of pullbacks of elements 
of Xh for all h, is included in a finite dimensional space. We can therefore find 
8' such that for all i € I, all 8 < 8', all h and all pullbacks Uh to Si of elements 
of Xh we have: 

\\uh - U h Rluh\\v>(\Ti\) < e\\uh\\L2(\Si\)- (75) 

This estimate is then transported back to M. D 

Commutator estimates. The operator HhRt does not commute with the 
exterior derivative but the following results provide some remedies. They show 
that the commutator: 

[H h Rtd] = ILXd - dlbX, (76) 

enjoys properties similar to I — HhR^- It is worth recalling here that Uh com- 
mutes with d, as do pullback operations. 

As before we first obtain a local stability result: 
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Lemma 4.5 For each S there is a constant C such that for all h, all T G 7^" 
and all u £ L 2 (st(T)) we have: 

||n fc [^,d]«|| La (r)<C||«|| La(Bt(T )). (77) 

- Proof: The following computation in the reference macroelement Si is useful: 

Rfdu-dRfu = (I - iPi)du + R s (iPtdu) - d((l - iPi)u) -dR 5 {iPiu),{78) 
= (I-R 5 )((d^)Au). (79) 

It shows that u is not differentiated. Remark also that (dipi) Am is equal to zero 
on a neighborhood of |Tj|. 

Next we remark that if we restrict attention to what happens on (a small 
neighborhood of) the simplex T G T£ ', the definition of R h can be taken to be: 

R h = r t n ■ ■ ■ r t 2 r t ± ) (80) 

where only the simplexes T m neighboring T are used (i.e. T m n T ^ 0). Then 
we can write: 

[Rl d] = J2 ^t n • • ■ i?T m+1 fe , d]i?L-. • • • R5 Tl ■ (81) 

m— 1 

In this sum the term where T m = T gives zero on a neighborhood of |T|, whereas 
the others constitute a small number (independent of T and h) of products each 
containing the regularization R^, composed (left and right) with terms which 
do not deregularize. 

From this the lemma follows. □ 

The estimates for Hh[R s h , d] then build up exactly as for I — 11/, R 5 h . 

Theorem 4.6 For each 5 and for each u G L 2 (M) we have: 

lim||n fc [i#,d]u|| L a(iio=0, (82) 

Proposition 4.7 For each Banach space X which can be considered as a sub- 
space o/L 2 (M) such that the canonical injection X — > L 2 (M) is compact we 
have, for each S: 

\im\\Tl h [Rid]\\ x ^ LHM) =0. (83) 

Proposition 4.8 For each e there is 5' such that for each S < 5' we have, for 
all h and all Uh G Xh : 

\\Il h [R s h ,d]u h \\ L 2 iM) <e|K|| L2(M) . (84) 

4.4 Applications 

While the properties of the operator 11/, R s h may seem weak - it does not com- 
mute with the exterior derivative, nor does it leave the Galerkin space invariant 

- we show that they are strong enough to prove some central estimates in numer- 
ical analysis. Specifically we are interested in estimates on discrete analogues 
of Hodge decompositions. 
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For any k the space of harmonic fc-forms is the space: 

n k {M ) = {a£ n k (M) : du = and W 6 f2 fc_1 (M) / u • dv = 0}. (85) 

On the other hand the space of discrete harmonic fc-forms is: 

Hh(M) = {u<=T k l (M) : dw = 0andVuer^- 1 (M) fu-dv = 0}. (86) 

For the rest of this section we consider a fixed (arbitrary) k, and put Xh = 
Y\(M). We will use also use the following notations: 

W h = {u h GX h : du h = 0} (87) 

V h = {«* 6 X h : \/w h eW h fu h -w h =0} (88) 

Similarly we let X be the completion of Vt k (M) with respect to the norm 
defined by: 

N&=Nlk + l|du|&. (89) 

Then we put: 

W = {uel : dw = 0} (90) 

K = {u e X ■. Vw e W fu-w = 0} (91) 

It will be useful to let Py denote the projection in X with range V and kernel 
W. 

The continuous Poincare-Friedrichs inequality is the assertion that there is 
a C > such that: 

VveV \\v\\v <C||du|| L 2. (92) 

Since the manifold M is supposed to be compact (without boundary) V is a 
subspace of the Sobolev space H 1 (M). In particular V is compactly embedded 
in L 2 , which is a variant of Rellich compactness. A major point is that in 
general Vh is not a subspace of V. However we will prove analogues of the 
Poincare-Friedrich and Rellich properties for the family of spaces (Vh). 

The following is a discrete Poincare-Friedrichs inequality. 

Proposition 4.9 There is C > such that for all h we have: 

Vv h £V h \\v h \\ L 2<C\\dv h \\ L 2. (93) 

- Proof: Pick Vh € Vh- We remark that dPyVh = dvh- Let then: 

v h - Pvv h = du + f, (94) 

be the Hodge decomposition of Vh — PyVh (see e.g. Taylor [35] Chapter 5, §8, 
Proposition 8.2). That is, / is harmonic and the only thing we need to know 
about u is that its norm with respect to a space compactly embedded in L 2 , 
is controlled by the L 2 norm of Vh ~ PvVh, hence also by the L 2 norm of Vh- 
Remark that dHhR^u + Hhf € Wh- In the following, all unspecified norms are 



L norms. We write: 






Ikll < K-dn fc i^«-n h /||, 




(95) 


< \\P v v h +du + f-dUhR S h u-n h f\\, 




(96) 


< \\P v v h + du- IL h R s h du + tt h R 5 h du - dtt h R s h u + f 


-Hfc/H, 


(97) 


< \\P v v h \\ + \\(I- n h R S h )du\\ + \\Uh[R s hl d]u\\ + ]](/- 


•n A )/||. 


(98) 
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We treat separately the three last terms on the right hand side. 

1) Concerning the second term we remark that, by the definition ([M]) : 

||(I-nfcitf)du|| < \\(I-IL h Ri)P v v h \\ + \\(I-Il h R 5 h )v h \\ + \\(I-Il h Ri)f\\. (99) 

In this estimate we have three terms. The first term is controlled by: 

\\(I-U h Ri)Pvv h \\ < \\(I -TLhRDWv^WPwhWx, (100) 



and we are ready to apply Proposition 14.31 For the second term we choose 
now 5 > such that the estimate of Proposition 14.41 holds with e = 1/2. The 
last term is controlled using Theorem 14.21 and the finite-dimensionality of the 
space of harmonic forms. Combining these three estimates we get, for a certain 
sequence e/, converging to 0: 

||(I - IIfcitf)d«|| < 1/2|M| + e h (\\v h \\ + \\dv h \\). (101) 

2) As already remarked, u is controlled in a space Y (e.g. the Sobolev space 
H 1 (M)) which is compactly embedded in L 2 . For the third term we can write: 

||n&[i^,d]«|| < \\IL h [Rid}\\ Y ^ v2 C\\v h \\ L 2, (102) 

and are ready to apply Proposition 14.71 

3) For the fourth term we use the smoothness of harmonic forms, and the 
finite-dimensionality of the space they form. 

Combining all three estimates we conclude that there is a sequence (e' h ) 
converging to such that: 

l/2|k|| < \\P v v h \\+e h (\\v h \\ + \\dv h \\) < (C + e' h )\\dv h \\+e' h \\v h \\. (103) 

This proves the proposition. □ 

An easy consequence which will be useful later is the following: 

Corollary 4.10 There is C > such that for all h we have: 

Vu h e X h inf \\uh - w/Jl 2 < C , ||du/ l || L 2. (104) 

w h €W h 

The discrete Poincare-Friedrich estimate also provides useful so-called Inf- 
Sup conditions for saddle-point problems, in the sense of Brezzi |10j . Consider 
for instance the following construction: A Fortin operator $/j : X — > Xh can be 
constructed by assigning to u G X the unique clement u% € Xh, solution of: 

u h eli / W h e X h Ju h -u h + J Ph ■ du' h = Ju-u' h , 1Q ^ 



p h e dX h \ Vp' h e dX h J p' h ■ du h = J p' h ■ du 

The announced Inf-Sup condition is: 

Corollary 4.11 There is C > such that for all h we have: 

■ , J'Ph -du h 

ml sup -^- — n — - — > 1/C. (106) 

Ph€dX h eXh \\p h L 2 \\Uh\\V> 
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This guarantees the stability (uniform boundedness) of the above Fortin oper- 
ators $h : X — > X. Notice that if u G H k (M) then d$;,.u = and by density: 



Vv G fifc (M) u-dv = 0, (107) 

so that $/j maps H k (M) into Hf^M). The Fortin operator can therefore be 
used to obtain estimates on discrete harmonic forms. We will use the following. 
We have: 

\\(I-$ h )\ HHM) \\x->x^0. (108) 

In particular $h : T-L (M) — > 7^(M) is eventually injective. Using the coho- 
mological properties of the preceding section we know that 7i k (M) and H k (M) 
have the same dimension, so that $>h '■ H k (M) — > TL k (M) is invertible. Denote 
by ^>h : TL^M) — > 1-L k (M) its inverse. The following elementary estimate will 
be used: 

\\I-*h\\ < ||/-*fc||/(l-||I-* h ||)-^0. (109) 

The following is a strengthening of the discrete Poincare-Friedrichs estimate: 

Proposition 4.12 For each e there is hi such that for each h < h' we have: 
Vv h e V h \\v h - PyVhWw < e||dv h || L ». (no) 

- Proof: In the following all unspecified norms are L 2 (_M) norms. Remark that 
since d(vh — PvVh) = we have: 

dn h R s h (v h - PVv h ) = U h [d,R s h ](v h - P v vh). (Ill) 

Pick e > 0. From Propositions 14.71 14.81 and 14.91 there is 5 > and hi such that 
for h < hi: 

\\dR h R{{v h - P v v h )\\ < e/(2C)\\dv h \\, (112) 

where C is the constant appearing in Corollary 14.101 This corollary provides 
Wh G Wh such that: 

\\dU h R s h (v h - P v v h ) - w h \\ < e/2\\dv h \\. (113) 

Actually, using Proposition 14.41 we will suppose that 8 was also chosen such 
that for all h and all Uh G X^: 

||(/-ILXWII<e/2|KI|. (H4) 

We remark that wu is orthogonal both to Vh and PyVh- Now we write: 

\\v h -P v v h \\ 2 (115) 

= f(v h -P v v h )-(v h -P v v h ), (116) 

< J(v h - P v v h ) ■ (v h - P v v h - U h Rl(v h - P v v h ) + w h ) + (117) 

e/2||v fc -iVwh||||dt;fc||, (118) 

< J( Vh -P vVh ).((I-U h R s h ){v h -P v v h ))+ (119) 

e/2||v h -iV«/.||||dt»fc||, (120) 

< \\v h - P V v h \\(\\(I - U h Ri)v h \\ + \\(I- U h R d h )P v v h \\ + e/2||d^||ijl21) 
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This gives: 

Ik - Pvv h \\ < ||(I- Il h R s h )v h \\ + || (/ - TL h Ri)P v v h \\ + e/2||d%||(122) 
< e/2(\\v h \\ + \\dv h \\) + \\I-U h Ri\\ v ^ L 2\\P v v h \\ x . (123) 

Applying Proposition 14. 31 we get, for h small enough: 

IK-/V«/.||<e(||t;/ l || + ||dt;fc||). (124) 

By Proposition 14.91 this completes the proof. D 

The following corollary is often referred to as discrete compactness: 

Corollary 4.13 Suppose (uh) is a sequence of elements Uh G Xfc such that 
\\uh\\h 2 an d l|dw/i,||L 2 are bounded. Suppose furthermore that: 

Vh \fu' h G X k ~ l J u h ■ du' h = 0. (125) 

Then (out of every subsequence) one can extract a subsequence converging in 
L 2 . 

- Proof: Indeed Uh can be decomposed as Uh = Vh + Wh where Vh G Vh and 
Wh G Wh- In fact u>h is discrete harmonic. Now PyVh is bounded in V which 
is compactly embedded in L 2 , so (PyVh) is relatively compact in L 2 . Moreover 
^hWh is harmonic and is bounded in a finite-dimensional space. One then 
concludes using Proposition 14. 121 and estimate (|109[) respectively. □ 

A Complexes and the long exact sequence 

We now recall some definitions and facts from homological algebra, and refer 
the reader to Lang [57] and Gelfand-Manin [3D] for thorough expositions. In 
this paper by a complex we mean a sequence A' of vectorspaces equipped with 
linear operators d k : A k — >• A k+1 called differentials and satisfying, for each k, 
d k+1 o d k = 0. The cohomology group K k A* is (the vectorspace) defined by: 

R k A' = (kerd k : A k -> A k+1 )/(imd k - 1 : A k ~ l -»• A k ) (126) 

Most often the index k in d k is dropped and the complex is represented by a 
diagram: 

► A k - 1 -^A k ^ A k+1 -> ■ ■ ■ , (127) 

where it is implicit that the arrows represent instances of the differential d. 

If A* and B' are two complexes, a morphism of complexes f : A* —> B', 
is a sequence of linear operators f k : A k — > B k such that the following diagram 

commutes: 

A k _>. A k+i 

f k I / fc +! | (128) 

B k _^ B k+i 

A morphism of complexes f*:A'—t B* induces a sequence or linear maps 
K k f : R k A' -)■ H fc BV We have that H fc (/' ° g') = (H fc /*) o (R k g'). 
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Given two morphisms of complexes f',g':A'—> B*, a homotopy operator 
from /* to g* is a family h* of linear operators h k : A k -4 B k ~ 1 such that: 

g k -f k = h k+1 d k + d k - 1 h k . (129) 

If /* and g* are homotopic in the sense that a homotopy operator exists from 
one to the other, they induce the same maps H k A' — > H k B*. 

One of the basic tools used in this paper is the following: Suppose we are 
given three complexes A' , B* and C* , and morphisms of complexes /* : A' — > 
B* and g* : B* -4 C , providing for each k a short exact sequence: 

-4 A fe -4 S fc -4 C* -)• 0. (130) 

Then (see e.g. Lang [271 chapter XX, §2) one can construct a Zong exaci sequence 
linking the cohomology groups: 

• • • K k A' -4 R k B' -4 H fe C" -4 H fc+1 J 4* -4 H fe+1 B' • • • . (131) 

Here the first arrow is H fe /*, the second is H fc g* and the third one - usually 
denoted S k - is constructed by the snake lemma. 

Remark A.l In particular, if we are given a morphism of complexes g m : B* — > 
C consisting of surjections and define A k to be the kernel of g k , we remark that 
A* is a subcomplex of B* called the kernel complex, which has trivial cohomology 
if and only if g* induces isomorphisms in cohomology H k B m — > H k C m . 

Suppose that the above situation holds for complexes A', B' and C* (with 
/* and g') and also for A\, B' and C* (with /* and g'). Suppose furthermore 
that we have morphisms of complexes a' : A' -4 A\, j3* : B* — > B' and 
7* : C* — s- C* intertwining /* with /* and g\ with g*. Explicitly we have 
commuting diagrams: 

A* J " B k "" C'\ I) 



- 4 o 


/o fe 


B k 


9 k o 


c k 


a k 1 
A k 


A fe 


I3 k I 
B k 


at 


c k 






(132) 


n k+1 A' 


' k+1 a' 1 


R k+i A . 


(133) 



Then the following diagram commutes: 

U k A . ^ h E k B . n- g Rkc . d 

R k a' 1 K k fi' i H fe 7 * I 



B Composition with compact operators 

Lemma B.l Let (A n ) be a sequence of continuous operators from a Hilbert 
space Y to a Banach space Z , such that, for some operator A: 

Vu e X lim A n u = Au. (134) 

n— >oo 
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Then A is continuous and for any Banach space X and any compact operator 
K : X — >• Y we have: 

lim \\AK-A n K\\ x ^ z = 0. (135) 

n— >oo 

- Proof: Continuity of A follows from the uniform boundedness principle, and 
the norm convergence follows from the fact that K can be approximated in the 
L(X, Y) norm by finite rank operators, see e.g. Rudin [36] (Chapter 4, exercise 
13). □ 
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